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Second-order QCD corrections to the thrust distribution 
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We compute the next-to-next-to-leading order (NNLO) QCD corrections to the thrust distribution 
in electron-positron annihilation. The corrections turn out to be sizable, enhancing the previously 
known next-to-leading order prediction by about 15%. Inclusion of the NNLO corrections signif- 
icantly reduces the theoretical renormalisation scale uncertainty on the prediction of the thrust 
distribution. 
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Three-jet production cross sections and related event 
shape distributions in e^e~ annihilation processes are 
classical hadronic observables which can be measured 
very accurately and provide an ideal proving ground for 
testing our understanding of strong interactions. The de- 
viation from simple two-jet configurations is proportional 
to the strong coupling constant, so that by comparing the 
measured three-jet rate and related event shapes with 
the theoretical predictions, one can determine the strong 
coupling constant a^. 

The theoretical prediction is made within perturbative 
QCD, expanded to a finite order in the coupling constant. 
This truncation of the perturbative series induces a the- 
oretical uncertainty from omitting higher order terms. 
It can be quantified by the renormalisation scale depen- 
dence of the prediction, which is vanishing for an all-order 
prediction. The residual dependence on variations of the 
renormalisation scale is therefore an estimate of the the- 
oretical error. 

So far the threejet rate and related event shapes have 
been calculated [l|, 0] up to the next-to-leading order 
(NLO), improved by a resummation of leading and sub- 
leading infrared logarithms [1, 01 and by the inclusion of 
power corrections [5!]. 

QCD studies of event shape observables at LEP Q are 
based around the use of NLO parton-level event genera- 
tor programs [3] . It turns out that the current error on as 
from these observables [8] is dominated by the theoreti- 
cal uncertainty. Clearly, to improve the determination of 
as , the calculation of the NNLO corrections to these ob- 
servables becomes mandatory. We present here the first 
NNLO calculation of the thrust distribution, which is an 
event shape related to three-jet production. 

The thrust variable for a hadronic final state in e^e~ 
annihilation is defined as @ 



T — max 



LO 


7* 


qqg 


tree level 


NLO 


7* 


qqg 


one loop 




7* 




tree level 




7* 


qqqq 


tree level 


NNLO 


7* 


^ qqg 


two loop 




7* 


qqgg 


one loop 




7* 


qqqq 


one loop 




7* 


qqqqg 


tree level 




7* 


^ qqggg 


tree level 



where pi denotes the three- momentum of particle i, with 
the sum running over all particles. The unit vector n is 
varied to find the thrust direction ut which maximises 
the expression in parentheses on the right hand side. 



TABLE I: Partonic contributions to the thrust distribution in 
perturbative QCD. 



It can be seen that a two-particle final state has fixed 
T — 1, consequently the thrust distribution receives 
its first non-trivial contribution from three-particle final 
states, which, at order as, correspond to three-parton fi- 
nal states. Therefore, both theoretically and experimen- 
tally, the thrust distribution is closely related to three-jet 
production. 

Three-jet production at tree-level is induced by the 
decay of a virtual photon (or other neutral gauge boson) 
into a quark-antiquark-gluon final state. At higher or- 
ders, this process receives corrections from extra real or 
virtual particles. The individual partonic channels that 
contribute through to NNLO are shown in Table HI All of 
the tree-level and loop amplitudes associated with these 
channels are known in the literature [l^, [ll|, [l^, [13 ■ 

For a given partonic final state, thrust is computed ac- 
cording to the same definition as in the experiment, which 
is applied to partons instead of hadrons. At leading or- 
der, all three final state partons must be well separated 
from each other, to allow T to differ from the trivial two- 
parton limit. At NLO, up to four partons can be present 
in the final state, two of which can be clustered together, 
whereas at NNLO, the final state can consist of up to 
five partons, such that as many as three partons can be 
clustered together. The more partons in the final state, 
the better one expects the matching between theory and 
experiment to be. 
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The two-loop 7* — > qqg matrix elements were de- 
rived in [To| by reducing all relevant Feynman integrals 
to a small set of master integrals using integration-by- 
parts and Lorentz invariance identities, solved 
with the Laporta algorithm |T6j . The master inte- 
grals [13 were computed from their differential equa- 
tions [15] and expressed analytically in terms of one- and 
two-dimensional harmonic polylogarithms [isj . 

The one-loop four-parton matrix elements relevant 
here [l^l were originally derived in the context of NLO 
corrections to four-jet production and related event 
shapes [l^, [l^]- One of these four-jet parton-level event 
generator programs [l^l is the starting point for our cal- 
culation, since it already contains all relevant four-parton 
and five-parton matrix elements. 

The four-parton and five-parton contributions to three- 
jet-like final states at NNLO contain infrared real ra- 
diation singularities, which have to be extracted and 
combined with the infrared singularities |21| present in 
the virtual three-parton and four-parton contributions to 
yield a finite result. In our case, this is accomplished by 
introducing subtraction functions, which account for the 
infrared real radiation singularities, and are sufficiently 
simple to be integrated analytically. Schematically, this 
subtraction reads: 
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where dcr^^^Q denotes the real radiation subtraction 
term coinciding with the five-parton tree level cross 
in all singular limits 



section da§j^^Q 



VS,1 
NNLO 



22l |. Likewise, 



is the one-loop virtual subtraction term coincid- 



v,i 



da 

ing with the one-loop four-parton cross section dcr^^^^Q 
in all singular limits [23| . Finally, the two- loop correction 
to the three-parton cross section is denoted by da^'^j^p. 
With these, each line in the above equation is individu- 
ally infrared finite, and can be integrated numerically. 

Systematic methods to derive and integrate subtrac- 
tion terms were available in the literature only to 
NLO Jli, [H], with extension to NNLO in special 
cases [2^. In the context of this project, we fully de- 
veloped an NNLO subtraction formalism [23], based on 
the antenna subtraction method originally proposed at 
NLO [13, • The basic idea of the antenna subtraction 
approach is to construct the subtraction terms from an- 
tenna functions. Each antenna function encapsulates all 
singular limits due to the emission of one or two unre- 
solved partons between two colour-connected hard par- 
tons. This construction exploits the universal factorisa- 
tion of phase space and squared matrix elements in all un- 
resolved limits. The individual antenna functions are ob- 
tained by normalising three-parton and four-parton tree- 
level matrix elements and three-parton one-loop matrix 



elements to the corresponding two-parton tree-level ma- 
trix elements. Three different types of antenna functions 
are required, corresponding to the different pairs of hard 
partons forming the antenna: quark-antiquark, quark- 
gluon and gluon-giuon antenna functions. All these can 
be derived systematically from matrix elements (28j for 
physical processes. 

The factorisation of the final state phase space into 
antenna phase space and hard phase space requires a 
mapping of the antenna momenta onto reduced hard mo- 
menta. We use the mapping derived in (29j for the three- 
parton and four-parton antenna functions. To extract the 
infrared poles of the subtraction terms, the antenna func- 
tions must be integrated analytically over the appropriate 
antenna phase spaces, which is done by reduction [30] to 
known phase space master integrals [3l|. 

We tested the proper implementation of the subtrac- 
tion by generating trajectories of phase space points 
approaching a given single or double unresolved limit. 
Along these trajectories, we observe that the antenna 
subtraction terms converge towards the physical matrix 
elements, and that the cancellations among individual 
contributions to the subtraction terms take place as ex- 
pected. Moreover, we checked the correctness of the sub- 
traction by introducing a lower cut (slicing parameter) on 
the phase space variables, and observing that our results 
arc independent of this cut (provided it is chosen small 
enough). This behaviour indicates that the subtraction 
terms ensure that the contribution of potentially singular 
regions of the final state phase space does not contribute 
to the numerical integrals, but is accounted for analyt- 
ically. A detailed description of the calculation will be 
given elsewhere [s^ . 

The resulting numerical programme, EERAD3, yields 
the full kinematical information on a given multi-parton 
final state. It can thus be used to compute any infrared- 
safe observable related to three-particle final states at 
0{al). As a first application, we present results for 
the NNLO corrections to the thrust distribution here. 
Leaving aside numerically small pure singlet-corrections 
(which come from 7* — > ggg and related final states, and 
appear first at NNLO), the theoretical expression for the 
thrust distribution through to NNLO can be expressed by 
three dimensionless coefficients {A,B,C), which depend 
only on T and not on the scale of the process or on cou- 
pling constants and quark charges. These are obtained as 
coefficients of the thrust distribution normalised to the 



tree-level cross section ctq for e"* 



qq, evaluated for 



as — as{Q), where Q is the centre-of-mass energy of the 
process: 

I da dA /as\'^dB /as\^dC 

dT ~ V2^/ dT V2^/ dr^V2^/ dT ' 
The experimentally measured thrust distribution 

1 da 



Chad dT 

can then be expressed in terms of A, B, C by expanding 
Chad around a^. If the renormalisation scale of as is 
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FIG. 1: Thrust distribution at Q = Mz- 
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FIG. 2: Relative scale-uncertainty on thrust distribution at 
different orders in perturbation theory. 



chosen to be /i 7^ Q, additional terms proportional to 
powers of hi{fj?/Q^) appear, which are again expressed 
in terms of A, B, C and of the coefficients of the QCD 
/3-function. 

In the numerical evaluation, we use Mz = 
91.1876 GeV and a^iMz) = 0.1189 Q. Figure [J displays 
the perturbative expression for the thrust distribution at 
LO, NLO and NNLO, evaluated for ^ = Q = Mz- It 
can be seen that inclusion of the NNLO corrections en- 
hances the thrust distribution by around (15-20)% over 
the range 0.03 < (1 — T) < 0.33. Outside this range, one 
does not expect the perturbative fixed-order prediction 
to yield reliable results. For (1 — T) > 0.33, the leading- 
order prediction vanishes due to kinematical constraints 
from having only three partons in the final state; at NLO, 
(1 — T) > 0.42 is kinematically forbidden, and the spher- 
ical limit T — > 0.5 is reached only for infinitely many 
partons in the final state. For T ^ 0, the convergence of 
the perturbative series is spoilt by powers of logarithms 
ln(l — T) appearing in higher perturbative orders, thus 
necessitating an all-order resummation of these logarith- 
mic terms [3, 0], and a matching of fixed-order and re- 
summed predictions [33j . 

To estimate the theoretical error inherent to the per- 
turbative prediction, we vary the renormalisation scale in 
the interval /i £ [Mz /2;2 Mz]- The relative uncertainty 
at each order 



S = 



maxp((T(^)) 



2(T(/i Mz) 

is displayed in Figure [21 It can be clearly seen how the 



inclusion of higher-order corrections stabilises the pre- 
diction, and the uncertainty S is reduced by about 30% 
between NLO and NNLO. The increase in uncertainty 
above 1 — T — 0.33 is due to the vanishing of the leading 
order contribution; the perturbative fixed-order descrip- 
tion and its theoretical error become unreliable beyond 
this point. 

In Figure [21 we compare the theoretical NNLO predic- 
tion for the thrust distribution to experimental data from 
the ALEPH experiment [s^l ■ Similar measurements were 
carried out by all LEP experiments [sH] ■ The NNLO cor- 
rection is positive and for the same value of as, yields 
a prediction that is larger than at NLO. This indicates 
the need for an improved determination of as from event 
shape data, which takes the newly computed NNLO cor- 
rections into account. Work on this is ongoing. 

In this letter, we presented the first NNLO calcu- 
lation of event shapes related to three-jet production 
in e+e" annihilation. We developed a numerical pro- 
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FIG. 3: NNLO thrust distribution compared to experimental 
data from ALEPH H. 
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gramme which can compute any infrared-safe observable 
through to ©(ttg), which we apphed here to determine 
the NNLO corrections to the thrust distribution. These 
corrections are moderate, indicating the convergence of 
the perturbative expansion. Their inclusion results in 
a considerable reduction of the theoretical error on the 
thrust distribution. Our results will allow a significantly 



improved determination of the strong coupling constant 
from jet observablcs. 
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